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Clear evidence for symplectic symmetry in low-lying states of 12 C and le O is reported. Eigenstates 
of 12 C and 16 0, determined within the framework of the no-core shell model using the JISP16 
AW realistic interaction, typically project at the 85-90% level onto a few of the most deformed 
symplectic basis states that span only a small fraction of the full model space. The results are nearly 
independent of whether the bare or renormalized effective interactions are used in the analysis. The 
outcome confirms Elliott's SU(3) model which underpins the symplectic scheme, and above all, points 
to the relevance of a symplectic no-core shell model that can reproduce experimental B(E2) values 
without effective charges as well as deformed spatial modes associated with clustering phenomena 
in nuclei. 



Recently developed realistic interactions, such as J- 
matrix inverse scattering potentials [l[ and modern two- 
and three-nucleon potentials derived from meson ex- 
change theory Q or by using chiral effective field theory 
Q , succeed in modeling the essence of the strong inter- 
action for the purpose of input into microscopic shell- 
model calculations that target reproducing characteris- 
tic features of light nuclei. The ab initio No-Core Shell 
Model (NCSM) [4| which employs such modern realistic 
interactions, yields a good description of the low- lying 
states in few-nucleon systems Q as well as in more com- 
plex nuclei like 12 C 0,0]. In addition to advancing our 
understanding of the propagation of the nucleon-nucleon 
force in nuclear matter and clustering phenomena 0, S] , 
modeling the structure of 12 C, 16 and similar nuclei 
is also important for gaining a better understanding of 
other physical processes such as parity-violating electron 
scattering from light nuclei @ and results gained through 
neutrino studies [lfj as well as for making better predic- 
tions for capture reaction rates that figure prominently, 
for example, in the burning of He in massive stars fill ] . 

In this letter we report on investigations that show that 
realistic eigenstates for low-lying states determined in 
NCSM calculations for light nuclei with the JISP16 realis- 
tic interaction [l[, predominantly project onto few of the 
most deformed Sp(3, R)-symmetric basis states that are 
free of spurious center-of-mass motion. This reflects the 
presence of an underlying symplectic sp(3, R) D su(3) D 
so (3) algebraic structure [221 ] . which is not a priori im- 
posed on the interaction and furthermore is found to re- 
main unaltered after a Lee-Suzuki similarity transforma- 
tion used to accommodate the truncation of the infinite 
Hilbert space by renormalization of the bare interaction. 
This in turn provides insight into the physics of a nucleon 
system and its geometry. Specifically, nuclear collective 



states with well-developed quadrupole and monopole vi- 
brational modes and rotational modes are described nat- 
urally by irreducible representations (irreps) of Sp(3,R). 

The present study points to the possibility of achieving 
convergence of higher-lying collective modes and reach- 
ing heavier nuclei by expanding the NCSM basis space 
beyond its current limits through Sp(3,R) basis states 
that span a dramatically smaller subspace of the full 
space. In this way, the symplectic no-core shell-model 
(Sp-NCSM) with realistic interactions and with a mixed 
Sp(3, R) irrep extension will allow one to account for even 
higher Ml configurations required to realize experimen- 
tally measured B(E2) values without an effective charge, 
and to accommodate highly deformed spatial configura- 
tions 12J that are required to reproduce a-cluster modes, 
which may be responsible for shaping, e.g., the second + 
state in 12 C and 16 O 0. 

We focus on the 0+, ground state and the lowest 
2 + (=2 ] F ) and 4 + (=4^) states in the oblate 12 C nucleus 
as well as the 0+ in the 'closed-shell' 16 nucleus. The 
NCSM eigenstates for these states are reasonably well 
converged in the N max = 6 (or 6hQ) model space with an 
effective interaction based on the JISP16 realistic interac- 
tion [H , which typically leads to rapid convergence in the 
NCSM evaluations, describes NN data to high accuracy 
and is consistent with, but not constrained by, meson ex- 
change theory, QCD or locality. In addition, calculated 
binding energies as well as other obscrvablcs for 12 C such 
as B(E2;2+^0+ ), B(Ml;l+^0+ ), ground-state proton 
rms radii and the l\ quadrupole moment all lie reason- 
ably close to the measured values. While symplectic al- 
gebraic approaches have achieved a very good reproduc- 
tion of low-lying energies and B(E2) values in light nuclei 
[HI, 14 1 and specifically in 12 C using phenomenological 
interactions [la ] or truncated symplectic basis with sim- 
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plistic (semi-) microscopic interactions 17 1, here, for 



the first time, we establish, the dominance of the sym- 
plectic Sp(3, R) symmetry in light nuclei, and hence their 
propensity towards development of collective motion, as 
unveiled through ab initio calculations of the NCSM type 
starting with realistic two-nucleon interactions. 

The symplectic shell model [l8|, is based on the 
noncompact symplectic sp(3,R) algebra. The classical 
realization of this symmetry underpins the dynamics of 
rotating bodies and has been used, for example, to de- 
scribe the rotation of deformed stars and galaxies [2p| . 
In its quantal realization it is known to underpin the 
successful Bohr-Mottclson collective model and has also 
been shown to be a multiple oscillator shell generaliza- 
tion of Elliott's SU(3) model. Consequently, symplectic 
basis states bring forward important information about 
nuclear shapes and deformation in terms of (A, /it), which 
serve to label the SU(3) irreps within a given Sp(3,R) 
irrep, for example, (0, 0), (A, 0) and (0, n) describe spher- 
ical, prolate and oblate shapes, respectively. 

The significance of the symplectic symmetry for a mi- 
croscopic description of a quantum many-body system 
emerges from the physical relevance of its 21 generators 
constructed as bilinear products of the momentum (p a ) 
and coordinate (qp) operators, e.g. p a pp, p a qp, and q a q/3 
with a, j3 = x, y, and z for the 3 spatial directions. Hence, 
the many-particle kinetic energy, the mass quadrupole 
moment operator, and the angular momentum are all 
elements of the Sp(3, M) D su(3) D so(3) algebraic struc- 
ture. It also includes monopole and quadrupole collective 
vibrations reaching beyond a single shell to higher-lying 
and core configurations, as well as vorticity degrees of 
freedom for a description of the continuum from irrota- 
tional to rigid rotor flows. Alternatively, the elements of 
the sp(3,R) algebra can be represented as bilinear prod- 
ucts in harmonic oscillator (HO) raising and lowering op- 
erators, which means the basis states of a Sp(3,R) irrep 
can be expanded in a 3-D HO (m-scheme) basis which 
is the same basis used in the NCSM, thereby facilitating 
calculations and symmetry identification. 

The basis states within a Sp(3,R) irrep are built by 
applying symplectic raising operators to a np-nh (n- 
particle-n-hole, n = 0, 2, 4, ...) lowest-weight Sp(3,R) 
state (symplectic bandhead), which is defined by the 
usual requirement that the symplectic lowering operator 
annihilates it. The raising operator induces a 2M7 lp-lh 
monopole or quadrupole excitation (one particle raised 
by two shells) together with a smaller 2ft£l 2p-2h correc- 
tion for eliminating the spurious center-of-mass motion. 
If one were to include all possible lowest-weight np-nh 
starting state configurations (n < N max ), and allowed all 
multiples thereof, one would span the full NCSM space. 

The lowest-lying eigenstates of 12 C and 16 were calcu- 
lated using the NCSM as implemented through the Many 
Fcrmion Dynamics (MFD) code 2l| with an effective 
interaction derived from the realistic JISP16 NN poten- 




FIG. 1: NCSM space dimension as a function of the maxi- 
mum TiQ. excitations, N max , compared to that of the Sp(3,R) 
subspace: (a) J = 0, 2, and 4 for 12 C, and (b) J = for 16 Q. 



tial [ij for different Ml oscillator strengths. For both 
nuclei we constructed all of the Op-Oh and 2MI 2p-2h 
(2 particles raised by one shell each) symplectic band- 
heads and generated their Sp(3, R) irreps up to N max = 6 
(6M7 model space). Analysis of overlaps of the symplec- 
tic states with the NCSM eigenstates for 2M7, 4M7, and 
QMl model spaces (N max = 2, 4, 6) reveals the dominance 
of the Op-Oh Sp(3,R) irreps. For the 0+ and the lowest 
2 + and 4 + states in 12 C there are nonncgligible overlaps 
for only 3 of the 13 Op-Oh Sp(3,R) irreps, namely, the 
leading (most deformed) representation specified by the 
shape deformation of its symplectic bandhead, (0 4), and 
carrying spin S = together with two (1 2) S = 1 ir- 
reps with different bandhead constructions for protons 
and neutrons. For the ground state of 16 O there is only 
one possible Op-Oh Sp(3,R) irrep, (0 0) S = 0. In addi- 
tion, among the 2M7 2p-2h Sp(3,R) irreps only a small 
fraction contributes significantly to the overlaps and it in- 
cludes the most deformed configurations that correspond 
to oblate shapes in 12 C and prolate ones in 16 0. 

The typical dimension of a symplectic irrep in the 
Nmax = 6 space is on the order of 10 2 as compared to 
10 7 for the full NCSM m-scheme basis space. As N max is 
increased the dimension of the J = 0, 2, and 4 symplec- 
tic space built on the Op-Oh Sp(3,R) irreps for 12 C grows 
very slowly compared to the NCSM space dimension (Fig. 
[TJl) . The dominance of only three irreps additionally re- 
duces the dimensionality of the symplectic model space, 
which remains a small fraction of the NCSM basis space 
even when the most dominant 2MI 2p-2h Sp(3,R) irreps 
are included. The space reduction is even more dramatic 
in the case of 16 O (Fig. QJ>). This means that a space 
spanned by a set of symplectic basis states is computa- 
tionally manageable even when high-M7 configurations 
are included. 

The overlaps of the most dominant symplectic states 
with investigated NCSM eigenstates for the 12 C and the 
16 in the 0, 2, 4 and 6M7 subspaces are given in TableU 
andHTl In order to speed up the calculations, we retained 
only the largest amplitudes of the NCSM states, those 
sufficient to account for at least 98% of the norm which 
is quoted also in the table. The results show that approx- 
imately 85% of the NCSM eigenstates for 12 C ( 16 0) fall 
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within a subspacc spanned by the few most significant 
Op-Oh and 2MI 2p-2h Sp(3,R) irreps, with the 2MI 2p- 
2h Sp(3, R) irreps accounting for 5% (10%) and with the 
leading irrep, (0 4) for 12 C and (0 0) for 16 O, carrying 
close to 70% (75%) of the NCSM wavefunction. 

Furthermore, the S = part of all three NCSM cigcn- 
states for 12 C is almost entirely projected (95%) onto only 
six S = symplectic irreps included in Table [H with as 
much as 90% of the spin-zero NCSM states accounted for 
solely by the leading (0 4) irrep. The 5 = 1 part is also 
remarkably well described by merely two Sp(3, R) irreps. 
Similar results are observed for the ground state of 16 O. 

Another striking property of the low-lying eigenstates 
is revealed when the spin projections of the converged 
NCSM states are examined. Specifically, as shown in 
Fig. [H their Sp(3, R) symmetry and hence the geometry 
of the nucleon system being described is nearly indepen- 
dent of the Ml oscillator strength. The symplectic sym- 
metry is present with equal strength in the spin parts of 
the NCSM wavefunctions for 12 C as well as 16 regard- 
less of whether the bare or the effective interactions are 
used. This suggests that the Lee-Suzuki transformation, 
which effectively compensates for the finite space trun- 
cation by renormalization of the bare interaction, does 
not affect the Sp(3,R) symmetry structure of the spa- 
tial wavefunctions. Hence, the symplectic structure de- 



TABLE I: Probability distribution of NCSM eigenstates for 
12 C across the dominant Op-Oh and 2MI 2p-2h Sp(3, R) irreps, 
HQ=15 MeV. 



0M1 2hQ 4MI 6MI Total 



J = 



S P (3, 



NCSM 



(0 4)5 = 46.26 12.58 4.76 1.24 64.84 

(1 2)5* = 1 4.80 2.02 0.92 0.38 8.12 

(1 2)5* = 1 4.72 1.99 0.91 0.37 7.99 

2hQ 2p-2h 3.46 1.02 0.35 4.83 

Total 55.78 20.05 7.61 2.34 85.78 

56.18 22.40 12.81 7.00 98.38 



Sp(3, 



NCSM 



(0 4)S = 46.80 12.41 4.55 1.19 64.95 

(1 2)5* = 1 4.84 1.77 0.78 0.30 7.69 

(1 2)5 = 1 4.69 1.72 0.76 0.30 7.47 

2hQ 2p-2h 3.28 1.04 0.38 4.70 

Total 56.33 19.18 7.13 2.17 84.81 

56.18 21.79 12.73 7.28 98.43 



Sp(3, 



NCSM 



(0 4)S = 51.45 12.11 4.18 1.04 68.78 

(1 2)5 = 1 3.04 0.95 0.40 0.15 4.54 

(1 2)5 = 1 3.01 0.94 0.39 0.15 4.49 

2hQ 2p-2h 3.23 1.16 0.39 4.78 

Total 57.50 17.23 6.13 1.73 82.59 

57.64 20.34 12.59 7.66 98.23 



tected in the present analysis for 6MI model space is what 
would emerge in NSCM evaluations with a sufficiently 
large model space to justify use of the bare interaction. 
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FIG. 2: Projection of the S = (blue, left) [and 5=1 (red, 
right)] Sp(3, R) irreps onto the corresponding significant spin 
components of the NSCM wavefunctions for (a) 0f 3 , (b) 2^, 
and (c) if in 12 C and (d) 0"j"„ in le O, for effective interaction 
for different HQ oscillator strengths and bare interaction. 



In addition, as one varies the oscillator strength Ml, 
the projection of the NCSM wavefunctions onto the sym- 
plectic subspace changes only slightly (see, e.g., Fig. [3]for 
the 0+ s state of 12 C and 16 0). The symplectic structure 
is preserved, only the Sp(3, R) irrep contributions change 
because the S = (S = 1) part of the NCSM eigen- 
states decrease (increase) towards higher Ml frequen- 
cies. Clearly, the largest contribution comes from the 
leading Sp(3,M) irrep (black diamonds), growing to 80% 
of the NCSM wavefunctions for the lowest Ml. These 
results can be interpreted as a strong confirmation of El- 
liott's SU(3) model since the projection of the NCSM 
states onto the i)Ml space [Fig. [3l blue (lowest) bars] is 
a projection of the NCSM results onto the SU(3) shell 
model. The outcome is consistent with what has been 
shown to be a dominance of the leading SU(3) symmetry 
for SU(3)-based shell-model studies with realistic inter- 
actions in 0MI model spaces. It seems the simplest of 



TABLE II: Probability distribution of the NCSM eigenstate 
for the J = ground state in 16 O across the Op-Oh and dom- 
inant 2hQ, 2p-2h Sp(3, R) irreps, ?ifi=15 MeV. 

mO. 2HQ, 6HQ, Total 



Sp(3,I 



NCSM 



(0 0)S = 
2MI 2p-2h 



50.53 



15.87 
5.99 



6.32 
2.52 



2.30 
1.32 



75.02 
9.83 



Total 



50.53 
50.53 



21.86 
22.58 



8.84 
14.91 



3.62 84.85 
10.81 98.83 
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FIG. 3: Ground + state probability distribution over 
OfiQ, (blue, lowest) to 6MI (green, highest) subspaces for 
the most dominant Op-Oh + 2htt 2p-2h Sp(3, R) irrep case 
(left) and NCSM (right) together with the leading irrep con- 
tribution (black diamonds), (0 4) for 12 C (a) and (0 0) for ie O 
(b), as a function of the hfl oscillator strength, N max — 6. 



Elliott's collective states can be regarded as a good first- 
order approximation in the presence of realistic interac- 
tions, whether the latter is restricted to a Qhil model 
space or richer multi-M7 NCSM model spaces. 

The 0+, and 2^ states in 12 C, constructed in terms 
of the three Sp(3, R) irreps with probability amplitudes 
defined by the overlaps with the NCSM wavefunctions 
for N max = 6 case, were also used to determine B(E2 : 
2+ 



0+ ) transition rates. The latter, increasing from 
101% to 107% of the corresponding NCSM numbers with 
increasing M7, clearly reproduce the NCSM results. 

In summary, we have shown that ab initio NCSM cal- 
culations with the JISP16 nucleon-nucleon interaction 
display a very clear symplectic structure, which is un- 
altered whether the bare or effective interactions for var- 
ious hil strengths are used. Specifically, NCSM wave- 
functions for the lowest 0+,, 2+ and 4+ states in 12 C and 
the ground state in 16 project at the 85-90% level onto 
a few Op-Oh and 2htt 2p-2h spurious center-of-mass free 
symplectic irreps. Furthermore, while the dimensionality 
of the latter is only w 10~ 3 % that of the NCSM space, 
they closely reproduce the NCSM B(E2) estimates. The 
wavefunctions for 12 C are strongly dominated by the 
three leading Op-Oh symplectic irreps, with a clear dom- 
inance of the most deformed (0 4)5' = collective con- 
figuration. The ground state of 16 is dominated by the 
single Op-Oh irrep (OO).!? = 0. The results confirm for 
the first time the validity of the Sp(3, R) approach when 
realistic interactions are invoked in a NCSM space. This 
demonstrates the importance of the Sp(3,R) symmetry 
in light nuclei while reaffirming the value of the simpler 
SU(3) model upon which it is based. The results further 
suggest that a Sp-NCSM extension of the NCSM may 
be a practical scheme for achieving convergence to mea- 
sured B(E2) values without the need for introducing an 
effective charge. In short, the NCSM with a modern re- 
alistic interaction supports the development of collective 
motion in nuclei which is realized through the Sp-NCSM 



and as is apparent in its 0MI Elliott model limit. 
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